Relation between the Liquid-Liquid Phase Transition 
and Dynamic Behavior in the Jagla Model 
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Using molecular dynamics simulations, we study a spherically-symmetric "two-scale" Jagla poten- 
tial with both repulsive and attractive ramps. This potential displays a liquid-liquid phase transition 
with a positively sloped coexistence line ending at critical point well above the equilibrium melting 
line. We study the dynamic behavior in the vicinity of this liquid-liquid critical point. We find that 
the dynamics in the more ordered high density phase (HDL) are much slower then the dynamics 
in the less ordered low density phase (LDL). Moreover, the behavior of the diffusion constant and 
relaxation time in the HDL phase follows approximately an Arrhenius law, while in the LDL phase 
the slope of the Arrhenius fit increases upon cooling. On the other hand, if we cool the system 
at constant pressure above the critical pressure behavior of the dynamics smoothly changes with 
temperature. It resembles the behavior of the LDL at high temperatures and resembles the behavior 
of the HDL at low temperatures. This dynamic crossover happens in the vicinity of the Widom line 
(the extension of the coexistence line into the one-phase region) which also has a positive slope. Our 
work suggests a possible general relation between a liquid-liquid phase transition and the change in 
dynamics. 



I. INTRODUCTION 



Water is the most important solvent for biological func- 
tion [l|, I2I , yet it possesses many properties that are not 
well understood. An open question of general interest 
concerning liquid water is the relation between a liquid- 
liquid (LL) phase transition |3,|j| and dynamic properties 
|a,|y,Qjla)llJjlll|- Based on analogies with other network 
forming hquids and with the thermodynamic properties 
of the amorphous forms of water, it has been suggested 
that, at ambient pressure, liquid water should show a 
crossover between fragile behavior at high T to strong 
behavior at low T H ll H HI [ll[ll| in the deep super- 
cooled region of the phase diagram below the homoge- 
neous nucleation line. This region, called the "no-man's 
land," may contain a LL critical point 3], the terminal 
point of a line of first order LL phase transitions. One 
current hypothesis concerns the possibility that water's 
anomalies are related to the existence of this LL critical 
point and its associated LL phase transition Hue [3, H 13 ■ 
Thus far there has been no direct experimental proof of 
this hypothesis, but recent experiments in nanoscale hy- 
drophilic pores [6j showing a line of sharp crossover in 
water diffusivity were interpreted in terms of the effects 
of the LL critical point. A dynamic crossover has also 
been associated with the LL phase transition in silicon 
and silica [ij, [ig . The relation between LL phase tran- 
sition and dynamics properties may be not limited to 
tetrahedral liquids, but is a general feature of all liquids 
near the LL critical point. Here we review the results of 
Ref. ^ on the Jagla model of liquid which models in- 
termolecular interactions using a spherically-symmetric 
soft-core potential. We also study phase segregation and 



dynamics below the LL critical point. 



II. TWO-SCALE SPHERICALLY-SYMMETRIC 
JAGLA RAMP POTENTIAL 

A spherically-symmetric potential with two different 
len gth scales has been studied recently [3, 113) Hit H^ ^^ 
El M m pi IM 111. Here, we study the hnear ramp 
potential with both attractive and repulsive parts [l9j. 
The potential is defined as 
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where Ur = 3.5C/o is the repulsive energy, Ua = ^ f^o is 

the attractive part, a is the hardcore diameter, b — 1.72a 

is the well minimum, and c = 3a is the cutoff at large 

distance [Fig. Q. 

We approximate the potential of Eq.(^ by a step func- 
tion with small steps AC/ = Uo/8 and implement the 
standard discrete molecular dynamics algorithm for par- 
ticles interacting with step potentials [2l|, |23, |23|, We 
use a as the unit of length, particle mass m as the unit 
of mass and Uq as the unit of energy. The simulation 
time is therefore measured in units of ay^m/Uo, temper- 
ature in unit of Uo/kB, pressure in units of Uo/a^ and 
density p in units of Na^/L^, where L is the size of the 
system and N — 1728 is the number of particles. We 
use the Berendsen thermostat and barostat for constant 
temperature and constant pressure simulations [29j. 
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FIG. 1: The "two-scale" Jagla ramp potential with attractive 
and repulsive ramps. Here Ur = 3.5f/o, Ua ~ —Uo, a is the 
hard core diameter, b — 1.72a is the soft core diameter, and 
c = 3a is the long distance cutoff. In the simulation, we use 
a as the unit of length, and Uo as the unit of energy. 
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FIG. 2: The equation of state P{T, p) for the ramp potential 
in the vicinity of the LL phase transition. Lines indicate P(T) 
isochores of 21 densities (L — 15.0, 15.2, 15.4, ..., 19.0 from top 
to bottom). The LL critical point (closed circle) is located at 
P — 0.243, and T = 0.375, corresponding to the point at the 
highest temperature isochore crossing. 



III. RESULTS 

The equation of state of the Jagla model [Fig. |5] is 
obtained by using constant volume and constant tem- 
perature simulations (NVT-ensemble) . The model dis- 
plays a LL critical point located at Pc = 0.243 ± 0.003, 
Tc = 0.375 ± 0.002, and pc = 0.37 ± 0.01. We determine 
the LL coexistence line by Maxwell construction on the 
isotherms [Fig. |3]. The coexistence line has a positive 
slope of 0.96 ± 0.02kBa^'^- According to the Clapeyron 
equation 
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FIG. 3: The Van der Waals equation of state of the Jagla 
potential. The positively sloped LL coexistence line is ob- 
tained by the Maxwell rule construction of the van der Waals 
equation of state. 



the entropy in the HDL phase is lower than entropy in 
LDL phase. Thus the HDL phase is more ordered than 
the LDL phase, which is opposite to the LL transition 
found in simulations for water [j| and silicon jl5j. The 
position of the melting line is estimated as the temper- 
ature at which the solid - liquid first order phase tran- 
sition occurs upon gradually heating a system initially 
consisting of a crystal configuration obtained by sponta- 
neous crystallization at low temperature. The equilib- 
rium melting line is about 15% below the upper limit 
melting temperatures |30l l31j . Very recently, it has been 
suggested that by shrinking the attractive and repulsive 
ramp parameters of the Jagla model it is possible to 
change the sign of the slope of the LL coexistence line 
as well as to move the position of the critical point into 
the super-cooled region as in water 24]. The stability 
limit is defined by the HDL and low LDL spinodal lines 
at which 



{dP/dp)T = 0. 



(3) 



(2) 



We examine the LL phase segregation along each iso- 
chore as we cool the system down in the two-phase region. 
We did not observe any signs of the phase segregation at 
temperatures higher than the spinodal temperature for 
a given density. In Fig. [SJ the isochores terminate at 
points where the phase segregation is clearly visible. We 
show the structure of the N = 1728 system upon phase 
segregation, for L = 16.6a at T = 0.285. The liquid 
separates into two types of liquids - HDL and LDL - 
with a clear phase boundary. The density of each phase 
is consistent with the density at the coexistent line. We 
present a clear visualization of the phase boundary by 
displaying the bonds between particles within the range 
d < 1.2a. In the LDL phase, there are practically no 
such particles, while in the HDL phase any particle has 
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FIG. 4: Phase segregation for system with L = 16.6 at T = 
0.285. In the top panel, we present bonds between particles 
within a distance d < 1.2, indicating that the left side is 
the HDL phase and the right side corresponding to the LDL 
phase. The bottom panel is the histogram of densities along 
the horizontal direction, with a high density on the left and 
lower density on the right. 



at least one neighbor within this range. This is reflected 
in the high peak in the radial distribution function which 
appears in the HDL phase [Fig.0]. 
We also study the diffusivity D 



D — lini 



6i 



(4) 



where (. . .)t' denotes an average over all particles and 
over all t' along constant pressure paths in three different 
regions: (i) P — 0.250 path a above the critical point in 
one-phase region [Fig. |5], (ii) P = 0.225 path /? [Fig. [5], 
below the critical point in the LDL phase, and (iii) P ~ 
0.225 path 7, in the HDL phase [Fig.E]. 

Below the LL critical point along path 7, the diffusiv- 
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FIG. 5: A sketch of the P ~ T phase diagram for the two- 
scale Jagla model. Upon cooling along path /3, the liquid 
can remain in the LDL phase (non-Arrhenius) before crossing 
the LDL stability limit (LDL spinodal line). Upon heating 
along path 7, the liquid may remain in HDL phase up to the 
HDL spinodal line. Thus one does not expect any dramatic 
change in the dynamic behavior along the path /? and 7 in this 
case. Upon cooling at constant pressure above the critical 
point (path a), a dynamic crossover occurs from LDL- like 
behavior at high-temperature side to HDL-like behavior at 
low-temperature side upon crossing the Widom line j2|. 



ity D exhibits a temperature behavior different from the 
behavior in the LDL phase along path /3. The behavior 
of the diffusion constant and relaxation time in the HDL 
phase follows an approximately Arrhenius law, while in 
the LDL phase the slope of the Arrhenius fit increases 
upon cooling. As we cool along path /?, we can achieve 
very low temperatures without phase segregation so that 
we can measure the dynamic properties of the LDL in a 
metastable region below the coexistence line. 

On the other hand, by heating the HDL phase along 
path 7, we can reach the stability limit of the HDL de- 
noted by a HDL spinodal line. As we cross the HDL 
spinodal at constant pressure, the density of abruptly in- 
creases and the diffusivity of the system becomes equal 
to that of the LDL. In contrast, along path a in the one 
phase region above the LL critical point, D exhibits a dy- 
namic crossover from the LDL-like behavior on the high 
temperature side of the Widom line to HDL-like behav- 
ior on the low temperature side. This dynamic crossover 
is an indication of crossing the Widom line above the LL 
critical point [3. 

We also present results for the relaxation time r along 
different paths, where t is defined as the time when the 
intermediate scattering function decays to 1/e value for a 
certain q (the first peak of the static structure factor) j3^ . 
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FIG. 6: Dynamic behavior for Jagla potential. The T-dependence of diffusivity along constant pressure paths: (a) path /3 with 
P — 0.225 < Pc (open diamonds), and path 7 with P = 0.225 < Pc (closed circles). D in the HDL phase differs from its 
behavior in the LDL phase, (b) path a with 0.250 > Pc- A dynamic crossover occurs when the Widom line is crossed upon 
cooling I3. (c) Relaxation time r as a function of 1/T for P — 0.225. Similar to D, the relaxation time r in the HDL (path 7) 
differs from r in the LDL (path (3). (d) Relaxation time r along path a with P = 0.250. A dynamic crossover, as found in D, 
occurs as the system changes from "LDL-like" behavior to "HDL-like" behavior. 



Along path 7 with P — 0.225 [Fig.|HIc)], r differs from its 
behavior in LDL phase along path (3 [Fig. EJd)]. Above 
the LL critical point, along path a with P = 0.250, we 
find that r displays a crossover from LDL-like behavior 
on the high-temperature side to HDL-like behavior on the 
low temperature side, which is similar to the crossover 
found in D. 



IV. 



SUMMARY 



In summary, we study a simple spherically symmetric 
two-scale Jagla potential with both repulsive and attrac- 
tive parts. The system displays a LL phase transition 
well above the equilibrium melting line. The phase seg- 
regation shows direct evidence of the coexistence of two 



metastable states - HDL and LDL phase. The dynamic 
behaviors in the LDL and HDL are different from each 
other below the critical pressure, and change from LDL- 
like behavior to HDL-like behavior upon cooling above 
the critical pressure, suggesting the association of the 
changes in dynamic behavior with the LL critical point 
similar to the changes observed in water in nan-pores |g] . 
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